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Abstract 

The nonperturbative modification of tlie evolution equation for QCD gener- 
ating functional is applied to the analysis of the factorial moments of the intrajet 
particle distributions in QCD. Nonperturbative preasymptotic correction to the 
factorial moments scaling (F-scaling) is calculated. The results are compared 
with the available experimental data. The data show the possible presence of 
such a nonperturbative effect. 
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1 Introduction. 



Physics of jets provides one of the most important testing grounds for practical 
apphcations of the theory of strong interactions, QCD. Perturbative evolution of 
a quark-gluon jet is now well investigated (see, for example, [Q). In particular, 
one of the striking predictions of perturbative QCD is an angular ordering in the 
intrajet evolution due to the color coherence in gluon emission. Nonperturbative 
aspects of jet evolution are considerably less clear. From experience of applying 
the sum rules to the analysis of characteristics of the heavy resonances it 
is clear that when virtuality of radiated parton is of the order of 1 GeV it is 
necessary to take into account nonperturbative corrections due to background 
vacuum fields giving rise to QCD vacuum condensate. Technically this results 
in introducing additional soft interaction of hard parton currents with the soft 
medium parametrized by the condensate fields. A similar situation arises in 
describing an evolution of a jet in some nontrivial (e.g. nuclear) medium, so 
that inelastic collisions with it lead to dissipation of energy and momentum from 
the hard parton component. Although physically these situations are rather 
different, we expect the appropriate theoretical formalisms describing both cases 
to be rather similar. 

In both necessity of introducing new interaction vertices lead to the 

occurrence of new dimensional parameters. This in turn results in the specific 
scaling violations in the description of jet evolution, in which the perturbative in- 
teraction vertices are dimensionless (contain ratios of energies or virtualities, see 
[Q). Physically this could be interpreted as a beginning of nonperturbative pro- 
cess of the QCD string formation which takes energy from the perturbative com- 
ponent and pumps it into nonperturbative degrees of freedom. The first model 
describing the jet evolution taking into account the nonperturbative energy loss 
was proposed by I.M.Dremin 0. It was based on an analogy with the physics of 
high-energy electromagnetic showers in a medium, where scaling-invariant evolu- 
tion due to photon bremsstrahlung and electron-positron pair creation is accom- 
panied by the scaling-violating scattering due to inelastic interactions of shower 
particles with the atoms of the medium leading to their ionization [Q. Later 
the corresponding modified QCD evolution equations were analytically solved in 
0, where expressions for parton multiplicities in quark and gluon jets and the 
energy loss in the perturbative hard component were calculated. The discussion 
of this method, as well as of the relevant Monte-Carlo calculations 0, where 
nonperturbative component was described through the effective QCD lagrangian 
can be found in the review paper [^. Serious defect of Dremin's equations 
is that there the color interference in intrajet evolution, known to be of crucial 
significance for the description of its characteristics, is not taken into account. 
The problem we are discussing in this paper is to build up a formalism providing 
an analytical description of intrajet evolution and taking into account both color 
coherence and nonperturbative energy loss. 
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This program was started in the previous work of the authors p], where the 
modified jet evolution equations were solved in the leading logarithmic approx- 
imation (LLA) and the general evolution equations for QCD generating func- 
tional taking into account the nonperturbative dissipative corrections were writ- 
ten down. In |Q] we considered only the simplest characteristics of the intrajet 
evolution such as the nonperturbative dissipative corrections to the intrajet par- 
ticle multiplicities and distributions. In this paper we turn to a more detailed 
analysis of the jet particle content by calculating the nonperturbative dissipative 
corrections to the factorial moments characterizing it. 

Let us remind that perturbative QCD predicts an asymptotic independence 
of the moments of the intrajet particle distributions on energy (KNO scaling) 
[|T|. Closely related to the KNO scaling is a so-called factorial moments scaling 
(F-scaling) that holds in the high energy limit. A very clear review on this sub- 
ject can be found, e.g., in where one can also find the relevant references (see 
also a recent review |jlO[)- The perturbative preasymptotic corrections to this 
scaling are due to the running of the coupling constant and thus proportional 



to \Jas{Q'^) |]TT| (and therefore to l/VlnQ^). Below we shall restrict ourselves 
to the case of the frozen coupling constant, so the above-mentioned perturbative 
preasymptotic corrections are absent and the analysis is focused on the nonper- 
turbative preasymptotic F-scaling violation due to dissipative effects in the jet's 
evolution. 

Let us note, that a modification of the evolution equations for QCD generating 
functional similar to those in |§] and introducing an additional inactivation of 
perturbative evolution were recently described in fl. 



2 Generating functional for dissipative evolu- 
tion 

It is well known, that a formalism of generating functional is indispensable 
in analyzing the properties of particle distributions in QCD jets. In particular, 
it provides a unique possibility of analytically accounting for such delicate prop- 
erties of QCD jet evolution as color coherence and resulting angular ordering in 
the parton cascades. In the previous paper |Q we have proposed a nonpertur- 
bative generalization of the perturbative evolution equation for the QCD gen- 
erating functional describing the intrajet particle distributions. A new element 
is to introduce a dimensional coupling between the perturbative partonic and 
nonperturbative soft modes resulting in additional "dissipative" damping of the 
perturbative intrajet evolution. The modified evolution equation in MLLA (mod- 
ified leading logarithmic approximation) for the generating functional G{p,9\u) 

*We are grateful to I.M. Dremin for bringing our attention to this paper 



3 



reads 



^^dlno"^ = jdx^lK{x) [G{xp,e\u)G{{l-x)p,d\u)-G{p,d\u)]-[^ j9^^ 

p p 

where p is an energy of the parton decaying into the two ones having energies 
xp and (1 — x)p correspondingly, 9 is an angle between the off-spring partons 
and /3 is a dimensional coupling constant accounting for the coupling to the 
nonperturbative modes. In Eq. (|l|) 7q = INcO-sIt^ is a running coupling, [N^ = 3 
is a number of colors), K{x) = 1/a; — 1 + x(l — x)/2 is an MLLA kernel for gluon 
splitting. The initial condition for this equation has the form G(pmm, 6*1111111^^) = u, 
thus fixing the borderline for the possibility of perturbative evolution. From 
the requirement of probability conservation we get another constraint for the 
generating functional G{p, 6\1) = 1. The Eq. (0) is a further generalization of the 
one in [Q, where only the case of a = 1 was considered. 

In the absence of dissipative effects (i.e. for /3 = 0) the evolution described by 
Eq. ([l|) is characterized by a specific scaling . Namely, the generating functional 
G{p, 6\u) depends only on a certain combination of p and 6 having the form y = 
\n{p9/Qo), where Qq = Pmindmin is a minimal transverse momentum considered 
in terms of the perturbative evolution. This results, in particular, in the fact that 
both KNO scaling for intrajet multiparticle distributions and the (equivalent to 
it in the limit of high energy) F-scaling [^J hold in the limit of ?/ — > cxo (not only 
p —>■ 00). However, this is no longer true when the additional "dissipative" term 
is taken into account. The presence of the new dimensional coupling (3 leads, 
generally speaking, to the explicit violation of the perturbative scaling. 

Below we shall consider only the effects of the first order in a parameter 
e = (/3/p)", so we assume that the following decomposition takes place: 

G{p,e\u) = Goiy\u)-EG,iy\u) (2) 

so that the perturbative scaling for Go{p,6) and Gi{p,6) still holds. 

Substituting the decomposition (H) into the evolution equation Eq. (P, we 
get in the zeroth and first order in e 



dGo 
dy 



J dx 7o' K{x) [Go{y + Ina;) Go{y + ln(l - x)) - Go{y)] (3) 



dGi 
dy 



j dx 7o K{x) 



■ Gi(i/ + ln(l-x)) , 

GoUz + lnx + 

1 — a; r 



^ / , /. \\Giiy + Inx) ^ . . 
+ God/ + ln(l - - ^i(^) 



+ (4) 



dy 

with the initial conditions Go(2/ = 0) = m and Gi{y = 0) = 0. Let us note, that 
physically the evolution equation in the parton branching angle 9 Eq. (|I]) and the 



evolution equations in y Eqs.(^,^) have very different meaning. The former is 
really following the evolution of the intrajet distribution proceeding by emitting 
the gluons into decreasing opening angles. The latter ones are describing the 
evolution of our snapshot view of the cascade at a certain level. The initial 
condition here is formulated not at the origin of the parton tree, but at its very 
end. The condition G'o(y = 0) = u actually means, that a parton having a 
minimal transverse momentum can not evolve at all. Larger y mean that there is 
some phase space for degrading down to Qq. In this sense the evolution in y goes 
in the opposite direction to the more intuitive angular one. At the same time 
using y as an evolution variable is crucial in formulating a scale-invariant picture 
of the perturbative cascade, so it will be this evolution that will be studied in the 
following. Let us also note, that we shall restrict our consideration to the simple 
case of pure gluodynamics and frozen coupling constant 70 = const. 

Before proceeding to solving Eqs.(0,^ it is necessary to understand what 
are the integration limits in these equations. The most natural procedure is 
to demand the positivity of the arguments of the generating functions entering 
Eqs.(^^. This is equivalent to the requirement of stopping the evolution com- 
pletely at pt = Qo (i-e. no particle could be emitted with pt < Qq). Then the 
integration over the energy fraction inherited by the daughter parton is restricted 
to the interval 

i-e-y 

j dx . . . = j dx . . . (5) 

e-v 

In the high energy limit of large ?/, where we expect our basic perturbative treat- 
ment to hold, one can make a replacement 

1 

J dx . . . = J dx . . . , (6) 


provided that the corresponding integrals are convergent. 

In general from Eq.@ it is clear, that this constraint becomes increasingly 
important as we are moving towards the borderline of perturbative evolution , 
which corresponds to a starting point in the evolution equation in y. Thus the 
accuracy of the predictions following from the evolution equation Eq.(|I|) is to a 
certain extent limited by the necessity of fixing the initial conditions at the "least 
perturbative" end of the intrajet evolution. 

3 Perturbative Contribution to Moments 

For reader's convenience we briefly remind in this Section some important for- 
mulas in the solution [|I^ of the zeroth approximation Eq.(^ for the factorial 
moments of the intrajet particle distributions in the frozen coupling constant 
case 7o = const. Further details of this solution can be found in 
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Let us start with computing the mean multiphcity 

dGo{y\u) 



My) 



du 



(7) 



u=l 



From Eq.(^ and taking into account that G'o(y|l) = 1 we get a differential 
equation on mean multiphcity 

7o / dxK{x){nQ{y + ln(a;)) + no(?/ + ln(l - x)) - %(?/)} (8) 



dy 

with an obvious initial condition no(0) = 1. The solution is 

no{y) = e^^ (9) 
where 7 is determined by a transcendent equation 

7 = llB{l) (10) 

1 

B{a) = j dxK{x){x'' + {l-xY -I}, a ^ 0, -1, -2, ... (11) 


The function B{a) is plotted in Fig.l. The solution of Eq.(|TD|) determines 7 as a 
function of 7q (and therefore as a function of a^). The function 7(05) is plotted 
in Fig. 2. 

It is well known that in the high energy limit the KNO-scaling, character- 
izing the intrajet evolution, is equivalent to the factorial moments F-scaling , so 
that in the expansion 

Go{y\u) = l + f^^^^^^F^')n,{yr = l + ±^^^^^ y ^ 00 (12) 

q=l 1- q=l 1- 

where F^'^^ are the factorial moments of the intrajet particle distribution, these 
factorial moments do not depend on y. After substituting the above expansion 
into Eq. (^ one obtains a quadratic recurrent relation for the factorial moments 

= %(^,. E C^r^i°^i^i-U(7m, 7(g -m)), g > 1, (13) 

7y ^ 7o-Dl7yJ rn=l 

where 

R{a, b) = J dxK{x)x''{l - xf = B{a, b + 2) + ^B{a + 2,b + 2), 

B{a,b) being the Euler beta-function. The resulting values of the factorial mo- 
ments are easily obtained by numerically solving Eq. (p^ . In Fig. 4 the values of 
Fq for some as are shown. 
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4 Mean multiplicity 



In this section we calculate the corrections to the mean multiplicity in the first 
order in e = so let us take n = no — erii, where ni is a correction term that 
has to be calculated from the correction to the generating functional in Eq. (0) 



dGi{y\u) 

My) = 



du 



(14) 



From Eq.(D we get a differential equation on the correction to mean multiplicity 

5^ = 7o dxK{x) I + ^--^^ nM] + Te- 

lls) 

The above equation is a linear inhomogeneous one, so its solution is a sum of 
a general solution to a homogeneous equation and a particular solution to the 
inhomogeneous one. The solution of the homogeneous equation is proportional 
to e^^y , where 71 is a solution of 

71 = 7o^(7i - ") (16) 

In this case we have taken the integration limits to be (0, 1) , which here can be 
done without loosing accuracy. From Eqs.( [T0|JT6| ), and the fact that B{a) is a 



decreasing function of a, there follows an important inequality 

max(7, a) < 7i < a + 7. (17) 

for a > so that the solution to the homogeneous equation is dominating the 
high energy asymtotics y 00. The function 7 is plotted in Fig. 2 as a function 
of as and in Fig.3 as a function of a. In the following we shall need only the 
high energy asymptotics of the solution, which is provided by the solution of the 
homogeneous equation e^'^^: 




so the correction is known up to the prefactor which is determined from the initial 
conditions and can be found only from the full solution of Eq.(|l^). From the 



above equation and Eq.(|T^ we see that the relative magnitude of the dissipative 
correction is enhanced with respect to a naively expected product of the factors ~ 
e~°^ and e"'''^, the latter obviously being a typical 1/no preasymptotic correction. 

5 Factorial moments. 

In this section we will calculate the nonperturbative corrections to the factorial 
moments of the intrajet particle distribution. Let us remind (see, e.g., []T]] and 
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0) that the factorial moments Fq are defined as coefficients in the expansion of 
the generating functional in a power series in the mean multiplicity 

°° iv - lY 

G{p, e\u) = l + Y^ ^ r^Fqn", Fi = 1. (19) 

q=l ^• 

To the first order in the small parameter e we get from Eq. (^) 

oo P(0) oo ^f(0)„ _|_ F(l)r)„ 

= 1 + E \nl(u -If-eY. ' ' , ^ ' nt\u - 1)^ (20) 

g=l q=l 1- 

where eF^^^ is a dissipative correction to the g-th factorial moment and we have 
explicitly isolated the nonperturbative correction to the mean multiplicity cal- 
culated in the previous paragraph. From Eq. (^) we immediately read off the 
expression for the first order correction to the generating functional Gi. 

q=i y- q=i y- 

where ^q{y) = {qF^^^rii + F^^^nQ)nl~^ and the "initial conditions" Fi = Fi'^ = 
1 give $i(?/) = ni{y). A complete calculation of ^q{y) is possible but very 
cumbersome. To the accuracy we are interested in it will be sufficient to work 
with the leading contribution to ^q{y), greatly simplifying the calculations and 
making the results transparent. 
The full equation on ^q{y) reads 

K(y) = 7o / dx K{x) ^^%iy + In x) + ^^^^y %iy + ln(l - x)) - <l>,(y)| + 



Wo C^F^^e^^'^-^'^y I dxK{x) | '^>»(^ + l^^) (i_^)7fa-m) ^ ^22) 

(1 x)^ J 

The homogeneous part of Eg . ([22|) is the same as for in the equation for rii Eg. ([TED, 
so its solution has the same form Gqe"'^^ . The solution to the inhomogeneous 
eguation on ^q{y) is a sum of some exponents in y (this immediately follows from 
the linearity of Eg.(|2^) and the structure of $g/S with q' < q). Proceeding this 
way we get $g ~ $,^6^'^"™'-''''^ {m < q), where the sign ~ means that $q includes 
the terms from ^me^'^~"^^"'^ ■ All such components together with the solution of 
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the homogeneous part of the equation give ~ e*^''"'^'''''^^'''^^, e'^'^^ {q > m > !)• 
From Eg . ([T7|) it is clear that the leading term has the form e^''"^^'^^"*''^^^. Let us 
define 

^q{y) = e««-^)^+^i)^ = *g e^'^ where 7g = (g - 1)7 + 71. (23) 

Substituting this Ansatz to Eq.(p2D we have 
<7-i „ 

m=l 

+ 72^gy"rfxK(x){x^«-" + (l-x)^«-"-l} , q>l (24) 
\&i = C^j. From Eq.(P^ we immediately get a recurrent relation on 
2 <?— 1 

VI/, = ^ C^vl/^Ff „ [R{lm - a, 7(g - m)) + R{^{q - m), 7^ - a)] . 

7? - 70^(7? - a) ^1 

(25) 

Using Eqs.(^) and ( PB| ) one can extract FJ^^ from \&q: 



7gy 



= - gFf e^^i-^)^ g > 1 (26) 

(for g = 1 it is obvious that F^^^ = 0). Let us further introduce a convenient 
notation Eg by 

= - ^' (27) 

SO that in the zeroth and first orders in e = (^^^ we get 

G=l + f2 ^- ^^Qoi ^ _ (28) 

The function Eg is plotted in Fig.5. 

In order to compare our final result with the existing experimental data on 
the factorial moments it is convenient to introduce another function pg 

pg = ^ ~ ^'^^''^ = (^y vl/,e(^-"-^)^ (29) 

\QoJ 

From this formula we see that the convenience of using pg is that it directly mea- 
sures the nonperturbative contribution to the factorial moments. From Eq. ([25|) 
we see, that our formalism predicts that Pg is positive and does not depend on q. 
Some plots of pg are shown in Figs. 6 and 7. 
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Let us now turn to the analysis of the existing experimental data on factorial 
moments of the intrajet particle distributions. From Eq.(^) it is clear that our 
main prediction is that if we substitute into Eq.(p9D the experimentally measured 
factorial moments F^^^ and compute the quantity 

^e.v=^l_F^^^v/p(o)^=p^,^^ (30) 

the difference between the experimental and purely perturbative factorial mo- 
ments as measured by p^^^ should (a) be positive and (b) have no dependence on 
the moment's rank q. 

We have computed the function p'^^^ for the experimental data on intrajet 
factorial moments from [Q. We see, that p^^^ is indeed positive thus signalling 
the presence of the contributions beyond the standard perturbative ones in the 
experimental data. In our treatment it is tempting to assume, that this differ- 
ence is due to the influence of the soft nonperturbative modes on the process of 
intrajet's evolution. From Fig. 6 we see, however, that the data show a clear 
dependence of p^^^ on q. This varies somewhat with varying a, but always stays 
quite pronounced. This shows that our treatment of the nonperturbative effects is 
actually incomplete. However, this is something that had to be expected, because 
here the small parameter is actually 57 [|l^, so for large enough q the correction 



is already much larger than the "leading" perturbative term. Therefore with our 
accuracy we can meaningfully calculate only the corrections to the two first mo- 
ments, and here the independence of g as a starting approximation does not seem 
unreasonable. Let us also mention, that in order to get a complete understanding 
of the possible nonperturbative contributions, it is desirable to analyze the corre- 
sponding contributions brought in by conventional hadronization physics in the 
Monte-Carlo approach (actually the local parton-hadron duality we use seems 
to work quite well, see, e.g., the recent review [0]) and by taking into account 
the corrections originated by the running coupling. Finally, let us look at the 
sensitivity of our results to the choice of the value of the frozen perturbative cou- 
pling Us and the parameter a parametrizing the coupling to the nonperturbative 
modes. In Fig. 6 we show ratio pq for 3 different values of as and in Fig. 7 we 
show the dependence of the ratio Pq on the value of the parameter a fixing the 
functional form of the dissipative interaction. We see, that both dependencies 
are reasonably smooth. 



6 Conclusions. 

By using the previously developed formalism [^] we have computed a nonper- 
turbative preasymptotic contribution to the factorial moments of the intrajet 
particle distribution violating the perturbative scaling. The presence of the ad- 
ditional nonperturbative contribution is supported by the existing experimental 
data on the factorial moments of the intrajet particle distributions. To arrive 
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at a final conclusion on the magnitude of the nonperturbative contribution one 
also has to take into account the contributions from hadronization stage (con- 
ventionally parametrized by some standard MC scheme), as well as those due to 
the running of the coupling constant. 

Acknowledgments. We are grateful to I.M. Dremin and I.V. Andreev for 
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Russian Foundation for Basic Research under Grant N 96-02-16210. 
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Figure captions 

Fig.l The function B{a) and a graphical calculation of 7 and 71 (see Eqs.(|TO|, 
|1|)) for as = 0.22 (7^ = 0.42) and a = 1. 

Fig. 2 7o, 7 and 71 — « with respect to obtained from Eqs.( p!0| , p!6D . 

Fig. 3 7i dependence on a. See Eqs.(^, |l^) 

Fig.4 for three values of a^. See Eq.(p!3|) 

Fig. 5 Eg for different and a. See Eqs. (|27| , |28|) 



Fig. 6 Pq dependence on a^, a = 1. Based on 0PAL[|T4| data for gluon jets. 



Fig. 7 Pq dependence on a for as = 0.22 Based on 0PAL[|T4| data for gluon jets. 
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